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ABSTRACT 


We report the discovery of eight new Milky Way companions in ~ 1,800 deg 2 of op¬ 
tical imaging data collected during the first year of the Dark Energy Survey (DES). 
Each system is identified as a statistically significant over-density of individual stars 
consistent with the expected isochrone and luminosity function of an old and metal- 
poor stellar population. The objects span a wide range of absolute magnitudes (My 
from — 2.2mag to — 7.4mag), physical sizes (lOpc to 170pc), and heliocentric distances 
(30kpc to 330kpc). Based on the low surface brightnesses, large physical sizes, and/or 
large Galactocentric distances of these objects, several are likely to be new ultra-faint 
satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood- 
based algorithm to search for and characterize stellar over-densities, as well as identify 
stars with high satellite membership probabilities. We also present completeness esti¬ 
mates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric 
distances in the first-year DES data. 


1. Introduction 

Milky Way satellite galaxies provide a unique opportunity to study the low-luminosity thresh¬ 
old of galaxy formation and to better connect the baryonic component of galaxies with the dark 
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matter halos in which they reside. Prior to the Sloan Digital Sky Survey (SDSS), the faintest known 
galaxies had luminosities of ~ 10 5 Lq, and it was clear that the population of twelve “classical” Milky 
Way satellites was orders of magnitude smaller than would be naively expected in the cold dark 
matter paradigm (Klypin et al. 1999; Moore et al. 1999). Over the past decade, systematic searches 
of wide-held SDSS imaging have revealed fifteen additional arcminute-scale, resolved stellar over¬ 
densities (Willman et ah 2005a,b; Zucker et ah 2006a,b; Belokurov et ah 2006, 2007, 2008, 2009, 
2010; Grillmair 2006, 2009; Sakamoto & Hasegawa 2006; Irwin et ah 2007; Walsh et ah 2007) that 
have been either photometrically classified or spectroscopically confirmed as gravitationally-bound 
“ultra-faint” galaxies (Kleyna et ah 2005; Munoz et ah 2006; Martin et ah 2007; Simon & Geha 
2007; Geha et ah 2009; Koch et ah 2009; Belokurov et ah 2009; Walker et ah 2009; Carlin et ah 
2009; Aden et ah 2009; Willman et ah 2011; Koposov et ah 2011; Simon et ah 2011; Kirby et ah 
2013). These ultra-faint galaxies are the smallest, least luminous, least chemically enriched, and 
most dark matter dominated galaxies in the known Universe. 

Since all known ultra-faint Milky Way satellite galaxies were discovered in SDSS, the census 
of these objects is almost certainly incomplete due to the partial sky coverage ( 14,000 deg 2 ) and 

photometric magnitude limit (95% complete to r ~ 22 mag) of that survey. While only 27 Milky 
Way satellite galaxies are currently known, extrapolations of the luminosity function suggest that 
hundreds of luminous Milky Way satellites remain to be found in current and near-future wide-held 
optical imaging surveys (Tollerud et ah 2008; Hargis et ah 2014; He et ah 2015). 

The Dark Energy Survey (DES) is in the process of imaging 5,000 deg 2 of the southern Galactic 
cap in five photometric bands (Abbott et ah 2005; Diehl et ah 2014). The deep photometry of DES 
(r ~ 24 mag) will enable the detection of the faintest known satellite galaxies out to ~120kpc 
(compared to the SDSS limit of ~50kpc), and more luminous satellite galaxies out to the Milky 
Way virial radius (Rossetto et ah 2011). We have completed an initial search of the first year of 
DES data and report here on the eight most significant dwarf galaxy candidates discovered therein 
(Table 1). Since the physical nature of these candidates cannot be definitively determined with 
photometry alone, we refer to them by their discovery coordinates. If these candidates are later 
confirmed to be Local Group galaxies, they should be renamed after the constellation in which 
they reside: DES J0335.6—5403 (Reticulum II), DES J0344.3—4331 (Eridanus II), DESJ2251.2 
—5836 (Tucana II), DES J0255.4—5406 (Horologium I), DES J2108.8—5109 (Indus I), DES J0443.8 
—5017 (Pictor I), DES J2339.9—5424 (Phoenix II), and DES J0222.7—5217 (Eridanus III). If any 
are instead globular clusters, they would be known as DES 1 through N. After the completion of 
this work, we learned that DES J2108.8—5109 was previously identified by Kim et al. (2015) in 
data from the Stromlo Milky Way Satellite Survey and designated as a likely star cluster, Kim 2. 


2. Data Set 

DES is a wide-held optical imaging survey in the grizY bands performed with the Dark 
Energy Camera (DECam; Flaugher et al. 2010; Diehl et al. 2012; Flaugher et al. 2015). The 
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Fig. 1.— Locations of 27 known Milky Way satellite galaxies (blue; McConnachie 2012a) and eight 
DES dwarf galaxy candidates (red) in Galactic coordinates (Mollweide projection). The coordinate 
grid shows the equatorial coordinate system with solid lines for the equator and zero meridian. 
The gray scale indicates the logarithmic density of stars with r < 22 from SDSS and DES. The 
large contiguous region in the northern equatorial hemisphere shows the coverage of SDSS (Aim 
et al. 2014). The full DES footprint is outlined in red, and is now partially filled in by a region of 
~ 1,600 deg 2 near to the Magellanic Clouds and a region of ~ 200 deg 2 overlapping with the SDSS 
Stripe 82 held along the celestial equator. Both fields were observed during the first year of DES 
and that compose the Y1A1 data set. 
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DECam focal plane comprises 74 CCDs: 62 2kx4k CCDs dedicated to science imaging and 12 
2kx2k CCDs for guiding, focus, and alignment. DECam is installed at the prime focus of the 4- 
nreter Blanco telescope at Cerro Tololo Inter-American Observatory. In this configuration, DECam 
has a hexagonal 2.2-degree-wide field-of-view and a central pixel scale of 0.263 arcseconds. The full 
DES survey is scheduled for 525 nights distributed over five years. Here, we consider data collected 
between 15 August 2013 and 9 February 2014 during the first year of DES. 

The first internal annual release of DES data (YlAl) consists of ~ 12,000 science exposures 
processed by the DES data management (DESDM) infrastructure (Gruendl et ah, in preparation). 1 
Most of the YlAl footprint is covered by 2 to 4 overlapping exposures, or “tilings”, in each filter. 
Single exposures in a tiling are 90s in griz and 45s in Y. Here, we rely on the g- and r-band 
images for photometry, and use the i-band for star-galaxy separation. 

The DESDM image processing pipeline consists of image detrending, astrometric calibration, 
nightly photometric calibration, global calibration, image coaddition, and object catalog creation, 
as recently summarized in Balbinot et al. (2015). We refer to Sevilla et ah (2011), Desai et ah 
(2012), and Mohr et ah (2012) for a more detailed description of the DES single-epoch and coadd 
image processing. The SExtractor toolkit is used to create object catalogs from the processed 
and coadded images (Bertin 2011; Bertin & Arnouts 1996). The YlAl data release contains a 
catalog of ~ 131 million unique objects detected in the coadd imaging which are distributed over 
1,800deg 2 . This area includes ~ 200 deg 2 overlapping with the Stripe 82 region of SDSS, as well as a 
contiguous region of ~ 1,600 deg 2 overlapping the South Pole Telescope (SPT) footprint (Carlstrom 
et ah 2011). The DES imaging in the SPT region is unprecedented in depth. Figure 1 shows the 
coverage of YlAl in Galactic coordinates. 

We selected stars from the YlAl coadd object catalog based on the spread_model quan¬ 
tity output by SExtractor (Desai et ah 2012). To avoid issues arising from fitting the point- 
spread function (PSF) across variable-depth coadded images, we utilized the weighted-average 
(■ wavg ) of the spread-model measurements from the single-epoch exposures. Our stellar sample 
consists of well-measured objects with \wavgspread-model J| < 0.003, flags-{g, r, i} < 4, and 
mag err-auto-{g,r,i\ < 1. We also removed objects for which the mag-psf and mag-auto mea¬ 
surements differ by more than 0.5 mag because this is indicative of poor object characterization. 

We estimated the stellar completeness on a statistical basis over the full YlAl footprint by 
creating a test sample of high stellar purity using a color-based selection of r — i > 1.7. We then 
applied the morphology-based star selection criteria above that uses imaging in only a single band 
to evaluate the stellar completeness for the test sample. This approach is unbiased since the two star 
selection criteria are orthogonal, and benefits from a large statistical sample that is representative 
of the full YlAl data set. The stellar completeness was found to be >90% to g ~ 22 mag and 
falls to ~ 50% by g ~ 23 mag. We validated this completeness estimate with matched spectroscopic 


1 http://data.darkenergysurvey.org/aux/releasenotes/DESDMrelease.html 
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data in the overlap region between Y1A1 and SDSS Stripe 82. Based on studies with the DAOPHOT - ’ 
software package optimized for photometry in crowded stellar fields, we do not expect this stellar 
completeness to be reduced in the vicinity of DES satellite galaxy candidates relative to the Y1A1 
footprint at large. 

For point-like objects and a well-estimated PSF, the SExtractor mag-psf variables are ex¬ 
pected to give the best measurement of stellar fluxes. However, due to the aforementioned dif¬ 
ficulties with PSF estimation on deep coadded images, we chose instead to use the magjauto 
measurements. The magjauto measurements are found to give a less biased estimate of flux when 
compared to a stellar calibration sample from Pan-STARRS (Schlafly et al. 2012). Measured mag¬ 
nitudes are extinction corrected using the E(B — V ) dust maps from Schlegel et al. (1998). The 
relative calibration uncertainties are estimated via stellar-locus regression (Kelly et al. 2014) and 
are found to be ~2% across the survey footprint. Uncertaities in the offsets between the DES 
photometric system and the AB system are estimated to be ~ 1%. 


3. Search Methods 

Ultra-faint galaxies are discovered as arcminute-scale over-densities of individually resolved 
stars. The YlAl stellar object catalog is of such quality and depth that numerous stellar over¬ 
densities are immediately apparent from a visual scan. Several of these over-densities are not 
associated with any known star cluster, globular cluster, or satellite galaxy. To formalize the 
process of identifying new candidate satellite galaxies, we applied both (1) a simple spatial binning 
algorithm to facilitate inspection of the stellar density field, and (2) a matched-filter maximum- 
likelihood technique. These complementary approaches validated one another and the resultant list 
of candidates was vetted by both methods. 


3.1. Stellar Density Maps 

Several independent searches of the stellar density field were conducted. One approach in¬ 
volved direct visual inspection of coadded images. Other searches used binned stellar density maps 
constructed from the coadd object catalogs. As an example, we detail below how one of these maps 
was built and analyzed. 

We began by spatially binning the stellar catalog into equal-area pixels using the HEALPix 
scheme (Gorski et al. 2005). 2 3 We considered HEALPix pixel sizes of ~0?06 (nside = 1024) and 
~0?11 (nside = 512) to optimize sensitivity to satellites possessing different angular sizes. Since 


2 http://www.star.bris.ac.uk/~mbt/daophot/ 

3 http://healpix.sourceforge.net 



the stellar density is greatly enhanced in regions of the Y1A1 footprint near the LMC and Galactic 
plane, we further grouped the stars into larger regions of ~ 13 deg 2 (nside = 16) to estimate the local 
field density of stars. We corrected the effective solid angle of each pixel using the survey coverage, 
as estimated by mangle as part of DESDM processing (Swanson et al. 2008). 4 Several conspicuous 
stellar over-densities were immediately apparent after this simple spatial binning procedure. 

We increase our sensitivity to ultra-faint satellite galaxies by focusing our search on regions 
of color-magnitude space populated by old, low-metallicity stellar populations (Koposov et al. 
2008; Walsh et al. 2009). As a template, we used a PARSEC isochrone corresponding to a stellar 
population of age 12Gyr and metallicity Z = 0.0002 (Bressan et al. 2012). Sensitivity to satellites 
at varying distances was enhanced by considering 20 logarithmically spaced steps in heliocentric 
distance ranging from 20kpc to 400kpc (distance moduli 16.5 < M — m < 23.0). For each step 
in distance, all stars within 0.2 mag of the isochrone in magnitude-magnitude space were retained 
while those outside the isochrone template were discarded. We then created a significance map for 
each ~ 13 deg 2 region by computing the Poisson likelihood of finding the observed number of stars 
in each map pixel given a background level characterized by the local field density. 


3.2. Matched-filter Maximum-likelihood Method 

The simple approach described above is computationally efficient and easily generalizable. 
However, a more sensitive search can be performed by simultaneously modeling the spatial and 
photometric distributions of stars and incorporating detailed characteristics of the survey (vari¬ 
able depth, photometric uncertainty, etc.). One way to incorporate this information is through a 
maximum-likelihood analysis (Fisher 1925; Edwards 1972). Likelihood-based analyses have found 
broad applicability in studies of Milky Way satellites (e.g., Dolphin 2002; Martin et al. 2008a). 
Here we extend the maximum-likelihood approach to a wide-area search for Milky Way satellites. 
Similar strategies have been applied to create catalogs of galaxy clusters over wide-held optical 
surveys (e.g., Rykoff et al. 2014). 

Our maximum-likelihood search begins by assuming that the stellar catalog in a small patch 
of sky represents a Poisson realization of (1) a field contribution including Milky Way foreground 
stars, mis-classified background galaxies, and imaging artifacts, and (2) a putative satellite galaxy. 
The unbinned Poisson log-likelihood function is given by 

log£ = -/A + ^(l -Pi), (1) 

i 

where i indexes the objects in the stellar sample. The value pi can be interpreted as the probability 


4 http://space.mit.edu/-molly/mangle/ 
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that star i is a member of the satellite, and is computed as 

_ A m 

Pl A m + bi' 


( 2 ) 


Here, u represents the signal probability density function (PDF) for the satellite galaxy and is 
normalized to unity over the spatial and magnitude domain, S ; specifically, f :iU u d5 = 1. The 
corresponding background density function for the field population is denoted by b. 

We define the richness, A, to be a normalization parameter representing the total number of 
satellite member stars with mass >0.1 Mq. In Equation (1), / = f,udS represents the frac¬ 
tion of satellite member stars that are within the observable spatial and magnitude domain of the 
survey, and /A denotes the expected number of observable satellite member stars. 5 Maximiz¬ 
ing the likelihood with respect to the richness implies /A = p t . This condition makes clear 
that the satellite membership probability for each star in the catalog is a natural product of the 
maximum-likelihood approach. These membership probabilities can be used to prioritize targeting 
when planning spectroscopic follow-up observations. Figure 2 highlights the use of membership 
probabilities to visualize a low-surface-brightness satellite galaxy candidate. 


To characterize a candidate satellite galaxy, we explore the likelihood of the data, V, as a 
function of a set of input model parameters, 6. The signal PDF is assumed to be separable into 
two independent components, 


u(Vi | 6) = u s (V S)i | 0 S ) x u c (V Cti | Q c ). (3) 

The first component, u s , depends only on the spatial properties, while the second component, u c , 
depends only on the distribution in color-magnitude space. 

We modeled the spatial distribution of satellite member stars with an elliptical Plummer 
profile (Plummer 1911), following the elliptical coordinate prescription of Martin et al. (2008a). 
The Plummer profile is sufficient to describe the spatial distribution of stars in known ultra-faint 
galaxies (Munoz et al. 2012b). The spatial data for catalog object i consist of spatial coordinates, 
T> s l = {aj, <5*}, while the parameters of our elliptical Plummer profile are the centroid coordinates, 
half-light radius, ellipticity, and position angle, 6 S = {ao, do, rh, e, </>}• 

We modeled the color-magnitude component of the signal PDF with a set of representative 
isochrones for old, nretal-poor stellar populations, specifically by taking a grid of isochrones from 
Bressan et al. (2012) spanning 0.0001 < Z < 0.001 and 1 Gyr < r < 13.5 Gyr. Our spectral data for 
star i consist of the magnitude and magnitude error in each of two filters, V c ^ = {<&, <7 9j j, r*, oy ; j}, 
while the model parameters are composed of the distance modulus, age, and metallicity describing 
the isochrone, 6 C = {M—m, r, Z}. To calculate the spectral signal PDF, we weight the isochrone by 


'Mangle maps of the survey coverage are used in the calculation of the observable fraction at each position in the 


sky. 
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Table 1. Detection of new satellite galaxy candidates in DES Y1A1 


Name 

Q: 2000 

^2000 

m — M 

Map Sig 

TS Scan 

TS Fit 

r h 

e 


Epi 


(deg) 

(deg) 


( CT ) 



(deg) 


(deg) 


DES J0335.6-5403 (Ret II) 

53.92 

-54.05 

17.5 

24.6 

1466 

1713 

0 io"*" 0 ' 01 

u - lu _o.oi 

0 6+ 0 ' 1 
u.D_o 2 

72 ± 7 7 

338.1 

DES J0344.3-4331 (Erill) 

56.09 

-43.53 

22.6 

23.0 

322 

512 

0 03 - * -0 ' 01 
U.LM-0 01 

0 19+ 016 
u - iy _0.16 

90+ 30 
yu —30 

96.9 

DES J2251.2-5836 (TucII) 

343.06 

-58.57 

18.8 

6.4 

129 

167 

u - iZ _0.03 

- 

- 

114.9 

DES J0255.4-5406 (Horl) 

43.87 

-54.11 

19.7 

8.2 

55 

81 

0 04 +0 - 05 
uu ^_0.02 

- 

- 

30.6 

DES J2108.8-5109 (Ind I) 

317.20 

-51.16 

19.2 

5.5 

- 

75 

o oio+ 0 002 
u.u±u_ 0 002 

- 

- 

26.6 

DES J0443.8-5017 (Pic I) 

70.95 

-50.28 

20.5 

7.1 

- 

63 

0 02 - *” 0 07 
u - u ^_o.oi 

- 

- 

19.1 

DES J2339.9-5424 (Phell) 

354.99 

-54.41 

19.9 

5.1 

- 

61 

0 02 - *" 0 01 
u - u ^_0.01 

- 

- 

19.4 

DES J0222.7-5217 (Erilll) 

35.69 

-52.28 

19.9 

5.4 

- 

57 

0 007~*~° 005 

u.uu/ _ 0 003 

- 

- 

8.9 


Note. — Best-fit parameters from the maximum-likelihood fit assuming the composite isochrone described in Section 3.2. Uncertainties are 
calculated from the the highest density interval containing 90% of the posterior distribution. “Map Sig” refers to detection significance of the 
candidate from the stellar density map search method (Section 3.1). “TS Scan” refers to the significance (Equation 4) from the likelihood scan 
using a Plummer model spatial kernel with half-light radius r ^ = 0? 1 (Section 3.2). “TS Fit” denotes the significance of the likelihood method 
using the set of best-fit parameters. Ellipticities and position angles are not quoted for lower significance candidates where they are not well 
constrained by the data. For objects with significant ellipticity, the half-light radius is measured along the elliptical semi-major axis. Hpi is the 
estimated number of satellite member stars with g < 23 in the stellar catalog. 


-54.0 


-54.1 



-54.2 


54.0 53.8 

<*2000 (deg) 


54.0 53.8 

02000 (deg) 


Fig. 2.— Left: False color gri coadd image of the 0?3 x 0?3 region centered on DESJ0335.6 
—5403. Right: Stars in the same field of view with membership probability pi > 0.01 are marked 
with colored circles. In this color map, red signifies high-confidence association with DES J0335.6 
—5403 and blue indicates lower membership probability. The membership probabilities have been 
evaluated using Equation (2) for the best-fit model parameters listed in Table 1. 
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a Chabrier (2001) initial mass function (IMF) and densely sample in magnitude-magnitude space. 
We then convolve the photometric measurement PDF of each star with the PDF of the weighted 
isochrone. The resulting distribution represents the predicted probability of finding a star at a 
given position in magnitude-magnitude space given a model of the stellar system. 

The background density function of the field population is empirically determined from a 
circular annulus surrounding each satellite candidate (0?5 < r < 2?0). The inner radius of the 
annulus is chosen to be sufficiently large that the stellar population of the candidate satellite does 
not bias the estimate of the field population. Stellar objects in the background annulus are binned 
in color-magnitude space using a cloud-in-cells algorithm and are weighted by the inverse solid angle 
of the annulus. The effective solid angle of the annulus is corrected to account for regions that are 
masked or fall below our imposed magnitude limit of g < 23 mag. The resulting two-dimensional 
histogram for the field population provides the number density of stellar objects as a function of 
observed color and magnitude (deg -2 mag -2 ). This empirical determination of the background 
density incorporates contamination from unresolved galaxies and imaging artifacts. 

The likelihood formalism above was applied to the Y1A1 data set via an automated anal¬ 
ysis pipeline. 6 For the search phase of the algorithm, we used a radially symmetric Plummer 
model with half-light radius rh = 0?1 as the spatial kernel, and a composite isochrone model con¬ 
sisting of four isochrones bracketing a range of ages, r = {12 Gyr, 13.5 Gyr}, and metallicities, 
Z = {0.0001,0.0002}, to bound a range of possible stellar populations. We then tested for a 
putative satellite galaxy at each location on a three-dimensional grid of sky position (0.7 arcrnin 
resolution; nside = 4096) and distance modulus (16 < M — m < 24; 16kpc to 630kpc). 

The statistical significance at each grid point can be expressed as a Test Statistic (TS) based 
on the likelihood ratio between a hypothesis that includes a satellite galaxy versus a field-only 
hypothesis: 

TS = 2 [log C(\ = A) - log £(A = 0)1 . (4) 

Here, A is the value of the stellar richness that maximizes the likelihood. In the asymptotic limit, 
the null-hypothesis distribution of the TS will follow a y 2 /2 distribution with one bounded degree 
of freedom (Chernoff 1954). We have verified that the output distribution of our implementation 
agrees well with the theoretical expectation by testing on simulations of the stellar field. In this 
case, the local statistical significance of a given stellar over-density, expressed in Gaussian standard 
deviations, is approximately the square root of the TS. 

6 The Ultra-faint Galaxy Likelihood (UGALI) code; detailed methodology and performance to be presented else¬ 
where. 
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4. Candidate Selection and Characterization 

The two search methods described in Section 3 each produce significance maps of the full 
Y1A1 footprint, where peaks in these maps represent the three-dimensional seed positions (a 2000 > 
^ 2000 j M — m ) of possible satellite galaxies. Seeds were selected from the union of the search 
methods. Statistical significance thresholds were set at > 5cr for the stellar density map method 
and TS > 45 for the matched-filter maximum-likelihood method, yielding ~ 50 seeds. Most of these 
were discarded as being attributed to steep gradients in the stellar density field, numerical effects 
near the survey boundaries, imaging artifacts, and large nearby galaxies resolved into multiple 
closely spaced catalog objects. For this reason, we did not pursue investigation at lower significance 
thresholds. 

The resulting seed list was compared against catalogs of known star clusters (Kharchenko 
et al. 2013; Harris 1996, 2010 edition) and Milky Way satellite galaxies (McConnachie 2012a) as 
well as catalogs of other astrophysical objects that can produce false positives, such as large nearby 
galaxies (Nilson 1973; Corwin 2004) or galaxy clusters (Rykoff et al. 2014). Associated seeds include 
the Reticulum globular cluster, the Phoenix dwarf galaxy, AM 1, NGC 1261, NGC 1291, NGC 1553, 
NGC1851, NGC 7089, NGC 7424, ESO 121-SC003, and ESO 201-SC 010. 

We explored the multi-dimensional parameter space for each unassociated seed using the emcee 
module for Markov Chain Monte Carlo (MCMC; Foreman-Mackey et al. 2013), 7 and the likelihood 
function described in Section 3.2 with flat priors on each of the input parameters. For each seed, 
we ran an MCMC chain with 100 walkers that each make 1000 steps including a burn-in stage of 
50 steps. This is sufficient to sample the region of parameter space near the maximum-likelihood 
estimate. Only seeds with well-constrained posterior distributions enter our candidate list of new 
Milky Way companions. 

Table 1 presents the eight most significant stellar over-densities in the Y1A1 data set consistent 
with being previously unknown dwarf galaxies. When comparing the significances obtained with 
the map-based and likelihood scan algorithms, it is worth noting that the two methods were applied 
assuming different size scales for the target satellites, and that kernel assumed for the likelihood 
scan (r/j = 0?1) is larger than the majority of candidates listed in Table 1. After fitting the spatial 
parameters of the candidates, all are detected with high significance using the likelihood-based 
method. The dependence of detection efficiency on assumed kernel size is discussed is Section 5.2. 

The physical characteristics of these objects, as determined by the follow-up MCMC likelihood 
analysis, are summarized in Table 2. The best-fit values and uncertainties are determined from the 
peak of the posterior distribution and the 90% highest posterior density interval (Box &; Tiao 1973). 
A significant correlation was observed between the age of the stellar isochrone and the heliocentric 
distance — a degeneracy that may be expected given the evolution of the main sequence turnoff. 


'emcee v2.1.0: http://dan.iel.fm/emcee/current/ 
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For some DES candidates, the posterior distribution for the distance is multi-modal. The distance 
estimates provided in Table 2 indicate the peaks in the posterior distribution. 

To compare with previously known Milky Way satellite galaxies, we convert from DES g- and 
r-band magnitudes to visual magnitudes using 

ODES = 9SDSS - 0.104(g , s D ss - esdss) + 0.01 mag 

'/"'DES = ?’SDSS - 0.102(g S DSS - ^sdss) + 0.02 mag (5) 

=> V = 5des - 0.487(5des - Does) - 0.025 mag. 

This transform from DES g and r magnitudes to E-band magnitudes was derived using an SDSS 
stellar calibration sample and the equations from Jester et al. (2005). The absolute magnitude 
of each satellite is calculated using the sampling formalism of Martin et al. (2008a). For bright 
satellites, this formalism yields a very similar estimate to the integration of the stellar luminos¬ 
ity function for the best-fit model. However, for fainter satellites, the uncertainty in the total 
magnitude can be dominated by shot noise arising from sparse sampling of the stellar population. 
In this case, the additional association of a single bright star can have a strong influence on the 
measured magnitude of a satellite. Similarly, the evolution of individual member stars can substan¬ 
tially change the total luminosity. To quantify the impact of shot noise on the derived luminosity 
estimates, we use a representative isochrone weighted by a Chabrier IMF to simulate an ensemble 
of satellites with similar characteristics to the observed candidates. The quoted uncertainty on 
the luminosity reflects the expected shot noise from stars in the magnitude range visible to DES, 
17 mag < g < 23 mag. 

The angular and physical half-light radii listed in Table 1 and Table 2 are both given as two- 
dimensional quantities. The deprojected (three-dimensional) half-light radius is a factor ~ 1.3 
larger than the projected half-light radius for a variety of common density profiles (Wolf et al. 
2010). For objects with measured ellipticity, we report the half-light radius measured along the 
semi-major axis. 

As illustrated in Figure 1, the DES candidates are distributed throughout the Y1A1 footprint 
and occupy a portion of the celestial sphere in the direction of the Magellanic Clouds where no ultra¬ 
faint galaxies were previously known. The DES candidates are widely distributed in heliocentric 
distance from ~30kpc (DES J0335.6—5403) to >300kpc (DES J0344.3—4331). 

5. Discussion 

Galaxies are distinguished from star clusters by having a dynamical mass that is substantially 
larger than the mass inferred from the luminous stellar population and/or a significant dispersion 
in metallicities indicative of multiple generations of star formation and a deep enough gravitational 
potential to retain supernova ejecta (Willman Sz Strader 2012). While devoted spectroscopic follow 
up observations are necessary to unambiguously classify these objects, the properties given in 
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Table 2 already provide strong clues as to which candidates are most likely to be galaxies. First, 
the large physical sizes of most of these objects are more consistent with the locus occupied by known 
satellite galaxies of the Local Group than with globular clusters of the Milky Way, as shown in 
Figure 3. All of the DES candidates are of comparable surface brightness to the ultra-faint galaxies 
detected in SDSS (McConnachie 2012a). The two most compact systems, DES J2108.8—5109 and 
DES J0222.7—5217, fall in between the known ultra-faint galaxies and the faintest Milky Way star 
clusters, e.g., Koposov 1 and Koposov 2 (Koposov et al. 2007; Paust et al. 2014), Segue 3 (Belokurov 
et ah 2010; Fadely et ah 2011; Ortolani et ah 2013), Munoz 1 (Munoz et ah 2012a), Balbinot 1 
(Balbinot et ah 2013), and Kim 1 (Kim & Jerjen 2015). For the most significant DES candidates, 
it is possible to estimate the ellipticity. Whereas globular clusters tend to have ellipticity < 0.2 
(van den Bergh 2008; Martin et ah 2008a), the best measured candidate, DES J0335.6—5403, has 
an ellipticity ~ 0.6, which is more consistent with the population of known ultra-faint galaxies. 

Further insight can be gained by fitting isochrones to the observed stellar distribution in 
color-magnitude space. Two independent maximum-likelihood implementations confirm that the 
DES candidates are generally consistent with old (r > 10 Gyr) and metal-poor stellar populations 
(Z 0.0002). The first of these analyses is the pipeline described in Section 3.2 used in a mode that 
varies age and metallicity in addition to spatial parameters and distance modulus in a simultaneous 
fit. The second color-magnitude fitting procedure adopts a similar likelihood formalism, but fits the 
spatial and photometric distributions of the stars in two separate phases. Instead of assuming an 
initial mass function, the second method weights the stars according to their proximity to the best- 
fit centroid location, and then evaluates the consistency between each star and a given isochrone in 
color-magnitude space given the photometric uncertainty for that star. The second method is more 
robust to complications that might arise from stellar incompleteness and/or imperfect modeling of 
the initial mass function. Age estimates and metallicity upper limits for four of the more significant 
DES candidates are reported in Table 2. Like the previously known ultra-faint dwarfs, the new 
DES systems are old and nretal-poor (e.g., Brown et al. 2014). The latter fitting procedure has 
also been applied to non-extinction corrected magnitudes to independently validate the extinction 
values from Schlegel et al. (1998). 


5.1. Review of Individual Candidates 

Brief comments on the individual galaxy candidates are provided below, and spatial maps and 
color-magnitude diagrams for each candidate are provided in Figures 4-11. The rightmost panels 
of these figures show the satellite membership probabilities of individual stars that are assigned by 
the likelihood fit using a single representative isochrone with r = 13.5 Gyr and Z = 0.0001. Stars 
with high membership probabilities contribute most to the statistical significance of each candidate. 
The constellation designation, should these candidates be confirmed as dwarf galaxies, is listed in 
parenthesis. 



15 - 



Fig. 3.— Local Group galaxies (McConnachie 2012a) and globular clusters (Harris 1996, 2010 
edition) occupy distinct regions in the plane of physical half-light radius (geometric mean of the 
major and minor axes) and absolute luminosity. The majority of DES satellite candidates (red 
dots) are more consistent with the locus of Local Group galaxies (empty blue shapes) than with 
the population of Galactic globular clusters (black crosses). Several of the faintest globular clusters 
and systems of ambiguous classification are indicated with x marks: Koposov 1 and Koposov 2 
(Koposov et al. 2007; Paust et al. 2014), Segue 3 (Belokurov et al. 2010; Fadely et al. 2011; Ortolani 
et al. 2013), Munoz 1 (Munoz et al. 2012a), Balbinot 1 (Balbinot et al. 2013), PSO J174.0675- 
10.8774 / Crater I (Laevens et al. 2014; Belokurov et al. 2014), and Kim 1 (Kim &; Jerjen 2015). 
Dashed lines indicate contours of constant surface brightness at fi = {25, 27.5, 30} mag arcsec^ 2 . 
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Fig. 4.— Stellar density and color-magnitude diagrams for DES J0335.6—5403. Top left : Spatial 
distribution of stars with g < 24 mag that are within 0.1 mag of the isochrone displayed in the 
lower panels. The field of view is 1?5 x 1?5 centered on the candidate and the stellar distribution 
has been smoothed with a Gaussian kernel with standard deviation 0?027. Top center: Radial 
distribution of stars with g — r < 1 mag and g < 24 mag. Top right: Spatial distribution of stars 
with high membership probabilities within a 0?5 x 0?5 field of view. Small gray points indicate 
stars with membership probability less than 5%. Bottom left: The color-magnitude distribution of 
stars within 0.1 deg of the centroid are indicated with individual points. The density of the field 
within a 1 deg annulus is represented by the background two-dimensional histogram in grayscale. 
The red curve shows a representative isochrone for a stellar population with r = 13.5 Gyr and 
Z = 0.0001 located at the best-fit distance modulus listed in the upper left panel. Bottom center: 
Binned significance diagram representing the Poisson probability of detecting the observed number 
of stars within the central 0.1 deg for each bin of the color-magnitude space given the local field 
density. Bottom right: Color-magnitude distribution of high membership probability stars. 





















DES J0335.6—5403 (Reticulum II, Figure 4): As the nearest and most significant candi¬ 
date, DES J0335.6—5403 is highly conspicuous in the Y1A1 stellar density maps, with ~ 300 
member stars brighter than g ~ 23 mag. In fact, an over-density of faint stars at this position 
is even visible in the much shallower Digitized Sky Survey images, although it was not de¬ 
tected by Whiting et al. (2007) and other searches of photographic material. Note that like 
the previously known ultra-faint dwarfs, DES J0335.6—5403 very likely contains several blue 
horizontal branch stars identified by the likelihood procedure, two of which are relatively far 
from the center of the object. Given its luminosity, radius, and ellipticity, DES J0335.6—5403 
is almost certainly a dwarf galaxy rather than a globular cluster. As illustrated in Figure 3, it 
is significantly more extended than any known faint globular cluster, and its elongated shape 
would also make it an extreme outlier from the Milky Way cluster population. Among known 
dwarfs, DES J0335.6—5403 appears quite comparable to Ursa Major II (Zucker et al. 2006a; 
Munoz et al. 2010). DES J0335.6—5403 is only ~ 23 kpc from the LMC, and measurements 
of its radial velocity and proper motion will provide strong clues as to whether it originated 
as a Milky Way satellite or fell into the Milky Way halo as part of a Magellanic group. 

DES J0344.3—4331 (Eridanus II, Figure 5): At a distance of >330kpc, DES J0344.3—4331 
is nearly a factor of three more distant than any known outer halo globular cluster, and its half- 
light radius of ~ 170 pc is inconsistent with the sizes of globular clusters. It is therefore very 
likely that this object is a new dwarf galaxy. The color-magnitude diagram of DES J0344.3 
—4331 closely resembles that of another distant Milky Way satellite, Canes Venatici I, with a 
well-populated horizontal branch covering a wide range of colors (Zucker et al. 2006b; Martin 
et al. 2008b; Okamoto et al. 2012). Its large distance places DES J0344.3—4331 in a very 
intriguing range of parameter space for studying the quenching and loss of gas in dwarf 
galaxies. As has been known for many years, dwarf galaxies within ~ 250 kpc of the Milky 
Way and M31 are almost exclusively early-type galaxies with no gas or recent star formation, 
while dwarfs beyond that limit often have irregular morphologies, contain gas, and/or are still 
forming stars (e.g., Einasto et al. 1974; Blitz & Robishaw 2000; Grcevich &; Putman 2009; 


Table 2. Properties of DES satellite galaxy candidates 


Name 

Distance 

M* 

M v 

r l/2 

icgioM 

z 


(kpc) 

(10 3 M 0 ) 

(mag) 

(pc) 

logiot G 7 r ) 


DES J0335.6-5403 (Ret II) 

32 

2 6 +0 - 2 
zo -0.2 

-3.6 ± 0.1 

551= 

10.08 ± 0.21 

< 0.0003 

DES J0344.3-4331 (Erill) 

330 

S3RI 

-7.4 ± 0.1 

1721” 

10.10 ± 0.23 

< 0.0006 

DES J2251.2-5836 (TucII) 

58 

sil 

-3.9 ± 0.2 

120~*” 30 

iZU -30 

- 

- 

DES J0255.4-5406 (Horl) 

87 

o 4+ 1 - 4 
z -0.7 

-3.5 ± 0.3 

00 I 30 

9.96 ± 0.21 

< 0.0005 

DES J2108.8-5109 (Indl) 

69 

n Q - ! -9 - 4 
°’ 8 -0.4 

-2.2 ± 0.5 

1211 

- 

- 

DES J0443.8-5017 (Pic I) 

126 

2.8±™ 

-3.7 ± 0.4 

43+ 153 
* -21 

10.00 ± 0.16 

< 0.0004 

DES J2339.9-5424 (Phell) 

95 

9 R+ 1 - 2 
2 - 8 -0.7 

-3.7 ± 0.4 

33t 2 ° 

- 

- 

DES J0222.7-5217 (Erilll) 

95 

0 9+ 0 - 9 
uy —0.7 

-2.4 ± 0.6 

ills 

- 

- 


Note. — Uncertainties are calculated from the the highest density interval containing 90% of the posterior 
distribution. Stellar masses are computed for a Chabrier initial mass function. 
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Fig. 5.— Analogous to Figure 4 but for DES J0344.3—4331. A large number of stars, including 
several probable horizontal branch members, are present at magnitudes fainter than the g < 23 mag 
threshold of our likelihood analysis. This threshold was set by the rapidly decreasing stellar com¬ 
pleteness at fainter magnitudes. However, it is likely that extending to fainter magnitudes would 
cause the best-fit distance modulus of DES J0344.3—4331 to increase. Better constraints on the 
properties of DES J0344.3—4331 require the stellar completeness to be robustly quantified in this 
regime. 
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Spekkens et al. 2014). The next most distant Milky Way dwarf galaxy, Leo T, is gas-rich 
and has hosted recent star formation (Ryan-Weber et al. 2008; de Jong et al. 2008); deeper 
optical and H I imaging to search for neutral gas and young stars in DES J0344.3—4331 
has the potential to provide new insight into how gas is stripped from low-mass dwarfs and 
reveal the minimum mass for maintaining star formation over many Gyr. A radial velocity 
measurement will also shed light on if DES J0344.3—4331 has already passed close to the 
Milky Way or whether it is infalling for the first time. Given its distance, it is unlikely to be 
associated with the Magellanic Clouds. Like DES J0335.6—5403, DES J0344.3—4331 is clearly 
detected in Digitized Sky Survey images dating as far back as 1976. 

• DES J2251.2—5836 (Tucana II, Figure 6): DES J2251.2—5836 is the third new satellite with 
a large enough size (120 pc) to be tentatively identified as a dwarf galaxy with DES multi¬ 
band photometry alone. It has a similar luminosity to DES J0335.6—5403, but is a much lower 
surface brightness system and is a factor of ~ 2 farther away. DES J2251.2—5836 is ~ 19 kpc 
from the LMC and ~ 37kpc from the SMC, making it a strong candidate for another member 
of the Magellanic group. In the surface density map in the upper left panel of Figure 6, the 
outer regions of DES J2251.2—5836 appear elongated and distorted. However, these features 
are likely a result of noise rather than real distortions (Martin et al. 2008a; Munoz et al. 
2010). The distribution of likely member stars in the upper right panel is much rounder. The 
high detection significance of this object demonstrates the power of the likelihood analysis to 
simultaneously combine spatial and color-magnitude information. 

• DES J0255.4—5406 (Horologium I, Figure 7): DES J0255.4—5406, at a distance of ~ 87kpc, 
has only a sparse population of RGB and HB stars visible in the DES photometry (a hint of 
the main sequence turnoff may be present at the detection limit). Given the small number of 
stars visible in the DES data, deeper imaging and spectroscopy will be needed to characterize 
the system more fully. Its Plummer radius of ~ 60 pc establishes it as a likely dwarf galaxy, 
twice as extended as the largest globular clusters with comparable luminosities. DES J0255.4 
—5406 (perhaps along with DES J0443.8—5017 and DES J2339.9—5424; see below) is signifi¬ 
cantly farther from the Milky Way than any previously known dwarf galaxy with My < —4, 
suggesting that tidal stripping may not be needed to explain the low luminosities of the 
faintest dwarfs. On its own, the Galactocentric distance is not necessarily a good indicator 
of the past importance of Galactic tides in shaping the photometric and spectroscopic prop¬ 
erties of satellites. The most important factor is the peri-Galacticon distance, which is not 
yet known for the new satellites. DES J0255.4—5406 is ~ 40 kpc away from the Magellanic 
Clouds and a factor ~ 2 closer to them than to the Milky Way, making it a potential Magel¬ 
lanic satellite. If it is (or was) associated with the Magellanic group, it is possible that tides 
from the LMC could have been relevant to its evolution. A measurement of the systemic 
velocity will help clarify whether DES J0255.4—5406 is currently near apocenter, infalling, or 
associated with the Magellanic system. 
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Fig. 6.— Analogous to Figure 4 but for DES J2251.2—5836. 
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Fig. 7.— Analogous to Figure 4 but for DES J0255.4—5406. 


Membership Probability Membership Probability Membership Probability Membership Probability 


















































21 - 


• DES J2108.8—5109 (Indus I, Figure 8): We identify DES J2108.8—5109 with Kim 2 (Kim 

et al. 2015). DES J2108.8—5109 is one of the faintest ( My ~ —2.2 mag) and most compact 
{ni 12 pc) of the satellites discussed here. The object is visible in the coadded DES images. 

If it is a globular cluster, as argued by Kim et al. (2015), DES J2108.8—5109 is fainter and 
more extended than most of the other outer halo clusters, such as AM 1, Eridanus, Pal 3, 
Pal 4, and Pal 14. DES J2108.8—5109 is ~37 kpc from the SMC, ~ 55 kpc from the LMC, 
and ~ 69 kpc from the much more massive Milky Way, so it is more likely a satellite of the 
Milky Way than of the Magellanic Clouds. 

• DES J0443.8—5017 (Pictoris I, Figure 9): DES J0443.8—5017 has a large enough radius 
to be a likely dwarf galaxy, but the uncertainty on the radius measurement is large enough 
to make it also consistent with the globular cluster population. DES J0443.8—5017 has a 
prominent blue horizontal branch and hints of an elliptical shape, but fewer member stars 
are detected in the DES data. More accurate measurements of size and shape from deeper 
imaging, and kinematics and chemical abundances from spectroscopy, will be required to 
determine the nature of this object. The large distance of DES J0443.8—5017 places it far 
enough behind the Magellanic Clouds that it is less likely to be a Magellanic satellite than 
many of the other new discoveries. 

• DES J2339.9—5424 (Phoenix II, Figure 10): DES J2339.9—5424 is quite similar to DES J0443.8 
—5017, but slightly smaller and closer. Again, we cannot draw firm conclusions on its nature 
without additional data, At ~ 43 kpc from the SMC and ~ 65 kpc from the LMC, it is unclear 
whether DES J2339.9—5424 could plausibly be a Magellanic satellite. 

• DES J0222.7—5217 (Eridanus III, Figure 11): DES J0222.7—5217 is the most compact of 
the newly discovered objects in both angular and physical units. Along with DESJ2108.8 
—5109, it is the most likely of the new discoveries to be a distant globular cluster rather 
than a dwarf galaxy. However, they could also be consistent with an extension of the dwarf 
galaxy locus to fainter magnitudes and smaller sizes. Even though it is one of the lowest 
luminosity system identified in the DES data so far, its compactness gives it a relatively high 
surface brightness, and like DES J2108.8—5109 and DES J0443.8—5017, it is clearly visible in 
coadded images. However, only a handful of likely member stars are resolved at the depth 
of the YlAl data, and significantly deeper imaging will be needed to better constrain its 
physical properties and stellar population. 


5.2. Detection Completeness 

Given that no additional ultra-faint Milky Way satellite galaxies have been confirmed outside 
of the SDSS DR7 footprint, despite the large areas of sky subsequently observed by SDSS and 
Pan-STARRS, it is interesting that multiple candidates have been found within the comparatively 
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Fig. 8.— Analogous to Figure 4 but for DES J2108.8—5109. 




1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


Fig. 9.— Analogous to Figure 4 but for DES J0443.8—5017. 
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Fig. 10.— Analogous to Figure 4 but for DES J2339.9—5424. 
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Fig. 11. — Analogous to Figure 4 but for DES J0222.7—5217. 
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Table 3. Expected detection efficiencies for Milky Way companions in DES Y1A1 


Name 

My 

Distance 

r h 

Efficiency 

Efficiency 

Efficiency 

Efficiency 



(kpc) 

(deg) 

(r ext = 0?029) 

(r ext = 0?057) 

( r ext = 0°1) 

( r ext = r h ) 

DES-J0335.6-5403 (Ret II) 

-3.6 

32 

0.100 

1.00 

1.00 

1.00 

1.00 

DES-J0344.3-4331 (Eri II) 

-7.4 

330 

0.030 

1.00 

1.00 

1.00 

1.00 

DES-J2251.2-5836 (Tuc II) 

-3.9 

58 

0.120 

0.25 

0.98 

1.00 

1.00 

DES-J0255.4-5406 (Hor I) 

-3.5 

87 

0.040 

0.62 

0.78 

0.55 

0.78 

DES-J2108.8-5109 (Ind I) 

-2.2 

69 

0.010 

0.96 

0.69 

0.18 

0.97 

DES-J0443.8-5017 (Pic I) 

-3.7 

126 

0.020 

0.92 

0.74 

0.30 

0.89 

DES-J2339.9-5424 (Phe II) 

-3.7 

95 

0.020 

1.00 

0.96 

0.74 

0.99 

DES-J0222.7-5217 (Eri III) 

-2.4 

95 

0.007 

0.24 

0.06 

0.00 

0.28 

Segue 1 

-1.5 

23 

0.073 

0.72 

0.99 

0.99 

0.99 

Ursa Major II 

-4.2 

32 

0.267 

0.06 

0.97 

1.00 

1.00 

Bootes II 

-2.7 

42 

0.070 

0.94 

1.00 

1.00 

1.00 

Segue 2 

-2.5 

35 

0.057 

1.00 

1.00 

1.00 

1.00 

Willman 1 

-2.7 

38 

0.038 

1.00 

1.00 

1.00 

1.00 

Coma Berenices 

-4.1 

44 

0.100 

1.00 

1.00 

1.00 

1.00 

Bootes III 

-5.8 

47 

1.666 

0.00 

0.00 

0.00 

0.96 

Bootes I 

-6.3 

66 

0.210 

1.00 

1.00 

1.00 

1.00 

Sextans 

-9.3 

86 

0.463 

1.00 

1.00 

1.00 

1.00 

Ursa Major I 

-5.5 

97 

0.188 

0.00 

0.30 

0.90 

0.98 

Hercules 

-6.6 

132 

0.143 

0.86 

1.00 

1.00 

1.00 

Leo IV 

-5.8 

154 

0.077 

1.00 

1.00 

1.00 

1.00 

Canes Venatici II 

-4.9 

160 

0.027 

1.00 

1.00 

1.00 

1.00 

Leo V 

-5.2 

178 

0.043 

1.00 

1.00 

1.00 

1.00 

Pisces II 

-5.0 

182 

0.018 

1.00 

1.00 

1.00 

1.00 

Canes Venatici I 

-8.6 

218 

0.148 

1.00 

1.00 

1.00 

1.00 


Note. — Detection efficiencies are calculated from many realizations of satellites with the properties (luminosity My , distance, 
Plummer profile angular half-light radius r^) of a given ultra-faint galaxy/candidate as they would have been observed in DES 
Y1A1. The simulated satellites are uniformly distributed throughout the SPT region of the Y1A1 footprint, excluding regions 
of high stellar density near to the LMC, i.e., ~ 1,600 deg^. The rightmost columns list the detection efficiencies for extraction 
regions of different radii, r e xt- Here, a detection constitutes > 5cr stellar excess with g < 23 within the extraction region given 
the local density of the stellar field, after selecting stars that are consistent with the isochrone of an old and metal-poor stellar 
population at the satellite distance (i.e., following the map-based detection algorithm described in Section 3.1). The extraction 
region radii are choosen to reflect size scales used in the map-based search (r e xt = {0?029, 0?057}), likelihood scan (r e xt = 0?1; 
Section 3.2), and matched to the size of the satellite (r e xt = r h )• Data for previously known satellites are taken from references 
compiled by McConnachie (2012b). 
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small area explored by DES thus far. 8 Without definite classifications, it is difficult to incorporate 
the DES candidates into constraints on the luminosity function of Milky Way satellite galaxies. 
However, it is still possible to quantify the sensitivity of the first-year DES search using simple 
semi-analytic estimates of the completeness. 

First, we calculated the probability that each new satellite could have been detected in the 
Y1A1 data. We began by generating a large number of realizations of each galaxy candidate dis¬ 
tributed uniformly over the Y1A1 footprint. The candidates were modeled using radially symmetric 
Plummer profiles and the realizations included shot noise due to the limited number of stars ex¬ 
pected to be in the observable magnitude range of DES. We then applied the simple map-based 
detection algorithm described in Section 3.1 to evaluate the detection efficiency. To be “detected”, 
the satellite must possess at least 10 stars brighter than our imposed magnitude limit (g < 23) and 
a large enough surface brightness to pass the visual search selection criteria. Specifically, we consid¬ 
ered extraction of varying sizes and computed the Poisson probability of detecting n sa t e mte + Afield 
stars when expecting ngeid stars based on the local field surface density. We tested extraction regions 
with sizes corresponding to the pixel areas in the map-based search algorithm (r ex t = {0?029,0?057}; 
Section 3.1) and the kernel size from the likelihood scan (r ex t = 0?1; Section 3.2), as well as an 
extraction radius set to the angular half-light radius of the simulated satellite. When computing 
the local field density, we selected only stars along the isochrone at the distance of the satellite 
with r = 12Gyr and Z = 0.0002 (see Section 3.1). Table 3 summarizes the expected detection 
efficiencies for the DES candidates when applying a 5 a statistical significance threshold, as in our 
seed selection procedure for the map-based search. The results show that all of the DES candidates 
would have been identified over a substantial fraction of the Y1A1 footprint with non-negligible 
probability, and for several candidates, near certainty. 

Table 3 also shows that the detection efficiency is sensitive to the size of the extraction region. 
Extended systems such as DES J2251.2—5836 are unlikely to found using the smallest extraction 
regions considered here, whereas the reverse is true for compact systems such as DES J0222.7—5217 
and DES J2108.8—5109. This size dependence accounts for the low significance of the two most 
compact candidates, DES J0222.7—5217 and DES J2108.8—5109, in the likelihood scan (Table 1). 
After allowing their spatial extensions to be fit, the detection significances of these candidates 
increase to a level well above the our imposed threshold. 

For comparison, Table 3 also provides detection efficiency estimates for previously known 
ultra-faint galaxies (assuming they were located in the SPT region of Y1A1 instead of their actual 
locations). We find that all of the SDSS ultra-faint galaxies, with the exception of the highly 
extended Bootes III, could have been readily detected in YlAl. We attribute these high detec¬ 
tion efficiencies to the deeper imaging of DES relative to SDSS and note that DES J2251.2—5836, 
DES J0255.4—5406, DES J2108.8—5109, and DES J0222.7—5217 have a substantially reduced de- 


8 The classification of PSO J174.0675-10.8774 / Crater I as a globular cluster or dwarf galaxy is currently ambiguous 
(Laevens et al. 2014; Belokurov et al. 2014). 
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tection probability when the magnitude limit is raised to r < 22 mag, comparable to the stellar 
completeness limit of SDSS. 

Our Y1A1 search sensitivity can be quantified in a more general way by considering an ensemble 
of satellites spanning a range of luminosities, physical sizes, and heliocentric distances. Figure 12 
presents the discovery potential of our Y1A1 search expressed as the detection efficiency with respect 
to these galaxy properties, estimated by the same method described above with r ext = rh- Nearby, 
luminous, and compact objects have a high probability of being significantly detected whereas 
objects that are more distant, faint, and extended are less likely to be found. The detection 
threshold in the plane of physical size and luminosity is nearly parallel to contours of constant 
surface brightness, and is weakly dependent on the distance, provided that a sufficient number of 
stars are detected. 

Since the Y1A1 search procedure described in Section 3 is a combination of the map-based and 
likelihood-based search techniques, the actual completeness of our search is likely slightly higher 
than estimated here. We expect the likelihood method to be more sensitive to extended low surface 
brightness systems because it combines spatial and color-magnitude information simultaneously. 
As the depth of DES imaging increases (2 to 4 tilings in Y1A1 compared to 10 tilings planned 
after 5 years) and more advanced techniques are applied to separate stars and galaxies at faint 
magnitudes (e.g., Soumagnac et al. 2013; Fadely et al. 2012), we anticipate that lower surface 
brightness satellites will become accessible. Our present study is optimized for the detection of 
relatively compact (ro < 0?2) and radially symmetric stellar over-densities. The search for extended 
low surface brightness features in the stellar distribution will be the focus of future work. 


5.3. Total Number and Spatial Distribution of Milky Way Satellite Galaxies 

The discovery of eight new dwarf galaxy candidates in ~ 1,600 deg 2 of Y1A1 not overlapping 
with SDSS Stripe 82 is consistent with expectations from the literature (Tollerud et al. 2008; 
Rossetto et al. 2011; He et al. 2015; Hargis et al. 2014). By empirically modeling the incompleteness 
of SDSS, Tollerud et al. (2008) predicted that 19 to 37 satellite galaxies could be found over the full 
DES footprint. More recent estimates based on high-resolution N-body simulations (Hargis et al. 
2014) and semi-analytic galaxy formation models that include baryonic physics (He et al. 2015) 
predict ~ 10 new detectable satellite galaxies in DES. Large uncertainties are associated with each 
of these estimates due to weak constraints on the luminosity function in the ultra-faint regime. 
Additionally, as noted in Section 5.1, some of the DES candidates may be globular clusters or 
may be associated with the Magellanic Clouds. In the latter case, it becomes more challenging to 
directly compare our results to the predictions above, which assume an isotropic distribution of 
Milky Way satellite galaxies. 

A number of studies, beginning with Lynden-Bell (1976), note that many Milky Way satel¬ 
lite galaxies appear to be distributed on the sky along a great circle, indicating a planar three- 
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Fig. 12.— Sensitivity of the our Y1A1 search expressed as the detection efficiency as a function of 
satellite luminosity, physical size, and heliocentric distance. Each panel corresponds to a different 
heliocentric distance. Contours of constant surface density are indicated with dashed black lines: 
(i = {25, 27.5, 30} mag arcsec -2 . 
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dimensional structure rather than an ellipsoidal or spherical distribution. This great circle has a 
polar orientation relative to the disk of the Milky Way. The discovery of most of the SDSS ultra¬ 
faint dwarfs in the north Galactic cap region increased the apparent significance of this alignment. 
However, since the primary region surveyed by SDSS is located in the direction of this so-called vast 
polar structure (Pawlowski et al. 2012), the true anisotropy of the Milky Way satellite population 
is not yet clear. The next generation of deep wide-held surveys should be able to address this issue 
with wider sky coverage. 

In this context, it is interesting to consider the locations of the eight new satellites reported 
here, which may increase the known Milky Way dwarf galaxy population by ~ 30%. The thickness 
of the vast polar structure defined by Pawlowski et al. (2012) is 29 kpc, and we find that the DES 
satellites have a dispersion of ~ 28 kpc from this plane. This result is perhaps not surprising given 
their proximity to the Magellanic Clouds, which played a large role in defining the original Lynden- 
Bell plane. In fact, the entire Y1A1 search area (with the exception of Stripe 82) is located quite 
close to the previously known plane of satellites. Thus, any satellite galaxies identified in this data 
set are necessarily close to the plane, and a selection of eight random positions within this area would 
likely have a similar dispersion relative to the polar structure. A more quantitative characterization 
of the distribution of Milky Way satellites awaits the completion of the DES survey, including areas 
farther away from the vast polar structure, as well as future results from Pan-STARRS. 


6. Conclusions 

We report on the discovery of eight new dwarf galaxy candidates associated with the Milky 
Way and/or Magellanic Clouds found in ~ 1,800 deg 2 of imaging data collected during the first year 
of DES. These satellites span a wide range of absolute magnitudes (—2.2 to —7.4mag), physical 
sizes (10 pc to 170 pc), and heliocentric distances (30kpc to 330kpc). The projected positions of 
the DES candidates are in close proximity to the Magellanic Clouds, and it is possible that some 
may be associated with the Magellanic system. 

The nature of these systems cannot be conclusively determined with photometry alone. How¬ 
ever, judging from their low surface brightnesses, ellipticities, and/or large distances, it is likely 
that several are new dwarf galaxies, in particular, DES J0335.6—5403, DES J0344.3—4331, and 
DES J2251.2—5836. If spectroscopically confirmed, the DES candidates may become the first ultra¬ 
faint galaxies identified outside the SDSS footprint, and would significantly increase the population 
of Local Group galaxies in the southern hemisphere. The proximity of DES J0335.6—5403, at 
~ 30 kpc, suggests that it may be an interesting target for indirect dark matter searches using 
gamma-ray telescopes (e.g., Ackermann et al. 2014). The implications of these candidate galaxies 
for indirect dark matter searches are discussed in a separate paper (Drlica-Wagner et al. 2015). 

The second year of the DES survey was completed on 15 February 2015. In addition to filling 
in regions of non-uniform coverage in the western portion of the Y1A1 footprint, the second season 
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expands the DES survey to encompass over 4,000 deg 2 . The sensitivity to ultra-faint satellite 
galaxies achieved with first-year DES data already exceeds that of SDSS (Section 5.2), and over 
the next five years, DES is expected to make an important contribution to our understanding of 
the Milky Way environment. 

During the preparation of this manuscript, we were sent an independent study by Koposov et 
al. using publicly released images from the first year of DES. Koposov et al. (2015) identify nine 
candidate satellites of the Milky Way and/or Magellanic Clouds, including eight that overlap with 
candidates presented here. The candidate that Koposov et al. refer to as Grus 1 is located outside 
the Y1A1 footprint in a region that was observed during the first year of DES with good image 
quality, but which did not have sufficient coverage in all bands to enter the coaddition stage of the 
standard DESDM pipeline. Therefore, the stars that comprise Grus I are not in the coadd object 
catalog that was used for this analysis. 

We note that we have not used the coordinates shared by Koposov et al. as seeds in our analysis, 
nor have we tuned our search algorithms based on knowledge of the candidates reported in their 
work. The search methods presented here yield significance maps of the entire Y1A1 footprint, 
and the reported detections are the most significant points in the maps unassociated with known 
objects. While our final choice of 5cr for the significance threshold for reportable galaxy candidates 
was made after our knowledge of the results from Koposov et al., the threshold was chosen to provide 
as much timely information to the astronomical community as possible with minimal likelihood of 
false positives, rather than for agreement with the Koposov et al. detections. We conclude that the 
independent discovery of these Milky Way companions by two separate teams using distinct object 
catalogs and search algorithms strengthens the case for follow-up by the astronomical community. 
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